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Abstract 

We study the contact process with stirring on In this process, particles occupy 
vertices of Z'^; each particle dies with rate 1 and generates a new particle at a randomly 
chosen neighboring vertex with rate A, provided the chosen vertex is empty. Additionally, 
particles move according to a symmetric exclusion process with rate N. For any d and 
A, there exists Ac such that, when the system starts from a single particle, particles go 
extinct when A < Ac and have a chance of being present for all times when A > Ac- Durrett 
and Neuhauser proved that Ac converges to 1 as A goes to infinity, and Konno, Katori 
and Berezin and Mytnik obtained dimension-dependent asymptotics for this convergence, 
which are sharp in dimensions 3 and higher. We obtain a lower bound which is new in 
dimension 2 and also gives the sharp asymptotics in dimensions 3 and higher. Our proof 
involves an estimate for two-type renewal processes which is of independent interest. 

Keywords: interacting particle systems, contact process, contact process with rapid stirring 


1 Introduction 

The contact process on with birth rate A > 0 and stirring rate A > 0 is the Markov process 
on {0,1}^'^ with generator C = defined, for any function / : {0,1}^'* M 

that only depends on finitely many coordinates and any ^ G {0,1}^ , by: 


(4xi/)(o = • E ^ 


{x,p}CZ^: 

X'-^y 



if 2: = X, 

C(x) 

if z = y, 


otherwise, 


( 1 . 1 ) 


(4oit/)(0= E (/r"°)-/(0) + ^- E E ’ 


^(x) = l 


2d 


^(x)=l y^x 


where = 


i if 2 : = x, 
^{z) otherwise. 


ie{o,i} (1.2) 


and we write x ~ y if |x — y|i = 1. We interpret each vertex x G Z"^ as a region of space which 
can either be empty (represented by state 0) or contain a particle (state 1). The dynamics 
can be described as follows: 

^Department of Mathematics, University of British Columbia, anna.levit@math.ubc. ca 
^Johann Bernoulli Institute, University of Groningen, d.rodrigues.valesin@rug.nl 


1 




• £exc describes the motion of particles, which is according to an exclusion dynamics with 
rate N. This means that for each pair x ~ y with ^(x) = 1 and ^{y) = 0, the particle at 
X jumps to y with rate N. For each pair x ~ y with ^(x) = ^(y) = 1, it will be useful to 
think that the particles at x and y exchange positions with rate N (though this amounts 
to no change in the conhguration 

• ^cont describes the birth and death of particles, which is according to a contact process 
with rate A, described as follows. Each particle dies with rate 1 and gives birth with 
rate A; when a birth occurs, a new particle is placed at a position chosen uniformly at 
random among the neighbors of the parent (if the chosen position is not empty, no new 
particle is born). 

This process is then a mixture of two well-studied classes of interacting particle systems: 
the exclusion process, introduced by Spitzer in [11], and the contact process, introduced by 
Harris in [5]. It was first considered in [3] by Durrett in Neuhauser, who were motivated by the 
earlier work [2] by DeMasi, Ferrari and Lebowitz. In [2], the authors added exclusion dynamics 
to particles on Glauber-type spin systems and showed that, if the rate of the exclusion is taken 
to infinity, the system converges to the solution of an associated reaction-diffusion equation. 

Our interest will be to consider the contact process with stirring from the point of view of 
its extinction-survival phase transition. To explain what we mean by this, let us first consider 
the basic contact process (that is, take the above definition with = 0). Assume = Il{o} 
(the indicator function of the origin) and consider the probability (which depends on d and A) 

P [for all t there exists x such that = 1] . (1.3) 

This is non-decreasing in A. We say that the process dies out if (1.3) is zero and survives if 
it is positive. The phase transition for the contact process is the statement that there exists 
Ac(Z'^) S (0,oo) such that the process survives if and only if A > Ac(Z'^) (for a proof of this, 
and of all the facts we state here about the basic contact process, see [10]). 

Since attempted births on occupied vertices (which we call “collisions”) produce no new 
particles, it is easy to show that the process given by the total number of particles, 
is stochastically dominated by a branching process with birth rate A and death rate 1. This 
comparison yields Ac(Z'^) > 1, and in fact it is known that 

< A,(Z'') < 2 , d>l. 

When stirring is introduced, one can consider Ac(Z'^,A^), still defined as the infimum of 
the values of A for which (1.3) is positive. In [3], Durrett and Neuhauser proved that 

lim Ac(Z'^,iV) = 1, d>l. (1.4) 

N—>-oo 

This result roughly means that, as we take N ^ oo - thus allowing particles to move more 
and more between each birth or death event - collisions have less effect and the critical rate 
approaches that of the associated branching process. 

Regarding the rate of convergence in (1.4), Konno showed in [7] that for every d there 
exist positive constants c^, Cd such that, for all > 1, 


Ac(Z'^,A^) - 1 G < 


W ’ N 


Cd log N Cd log N 
N ’ N 


Cd Cd 
7V1/3 J Ari/3 


if d > 3; 
if d = 2; 
if d = 1. 


(1.5) 
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For d > 3, in [6], Katori improved this to 

——j—— < liminf A^(Ac(Z'^,iV) - 1) < limsup A^(Ac(Z'^, iV) - 1) < (1-6) 

Zdyld — Ij N^oo N^oo ^d 

where G{-,-) is the Green function of discrete-time, simple random walk on that is, if 
/3i, /32,... are independent and uniformly distributed on { 2 ; : z ~ 0}, then 

G{x, y) = l{x= 2 y} + lE[#{n > 1 : = V - x}], 

where denotes the cardinality of a set. More recently, Berezin and Mytnik proved that, for 
d>3, the upper bound in (1.6) is sharp: 

lim N . iUZ^N) - 1) = (1.7) 

N^oo Id 

We now state our main results. For a measurable set ACM, |A| denotes the Lebesgue 
measure of A. 

Theorem 1.1 Assume d > 2 and let (A^, denote the positions of two particles under 

exclusion dynamics on with rate 1 (per edge) and Aq ~ Bq. Suppose {Xn)n>i is a sequence 
satisfying 

limsup—j-- < 1. 

iv^oo 3^-E[|{t<iV: Ai~.Ba|] 

Then, if N is large enough, the contact process with birth rate Xn and stirring rate N dies 
out. 


The following helps relate this theorem to the bounds given earlier: 


Proposition 1.2 Let As,Bs be as in Theorem 1.1. Let A'^ and B( be independent, continuous¬ 
time simple random walks on which jump from each vertex with rate 2d, and Ag B'o- If 
d > 3, then 


E [|{t < 00 : A^ ~ Bt}\] = E [\{t < 00 : A( ~ B',}\] = \ (1.8) 

If d = 2, then 

E[|{s<t:A,^ija|] ^ E[\{s<t:A[^Bi}\] ^ J_ 

t^oo log t t^oo log t 27r 


(1.8) already appeared in [1]. (1.9), which presents a higher technical challenge, is new. By 
putting together Proposition 1.2 and Theorem 1.1, we re-obtain the sharp bound of Berezin 
and Mytnik for d> 3 and obtain a new bound for d = 2: 


Corollary 1.3 If d>3, then 


If d = 2, then 


lim inf iV- (Ac(Z'^,iV) - 1) > 

N^oo 


G(0,0) - 1 
2d 


N 


lim inf 
N^cx) logN 


(Ae(Z^iV)-l) > —. 


Itt 


( 1 . 10 ) 
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Our bound for d = 2 is of the correct order given in (1.5). No matching upper bound has been 
proved, but since both cases d = 2 and d > 3 are treated by the unified statement of Theorem 
1.1, and since the resulting bound is sharp in the second case, it is reasonable to expect that 
(1.10) is also sharp. 

In our proof of Proposition 1.2, we study the processes (^^ — Bs)s>o and {A'^ — i?')s>o. 
The latter is exactly a random walk on and the former is the same random walk except 
that its jump rates at the set of neighbors of the origin are modified, and in particular it can 
never reach the origin. Hence, denoting by Aq the set of neighbors of 0, 

(l{y4s — Bg G Ao})s>o and (11{A'^ — G A/q U {0}})s>o, 


are two-type (0 and 1) renewal processes, and the amounts of time they spend at stage 0 have 
the same distribution. 

Motivated by this, we consider the following more general setting. Let and V be 

positive random variables and consider three independent sequences of independent random 
variables. 


\ ... ~ f/W, , C/f \ C/f ... ~ Ho, Hi, H2, ... ~ I". 

For i G {1,2}, define sequences (Sn^)n>o as follows: 


Also define, for i G {1, 2} and t > 0, 


(i) 

Kf = 




(d c(d 1 
2n-|-lJ 


\n=0 


= t — 


|o.‘ln ( Uld2+i.s 


dd 1 

■’2n-|-2J 


\n=0 


( 1 . 11 ) 


= 0 and, for all n > 0, and " S^+i = K- (1.12) 


(1.13) 


We think of the sequence {Sn^)n>o as describing a renewal process which alternates between a 
“u^*^-state” (where it stays for an amount of time distributed as H^*^) and a “u-state” (where 
it stays for an amount of time distributed as V). Then, represents the total amount of 
time spent in the u^-^^-state before instant t. We prove: 


Theorem 1.4 E[(C/(2))2] < oo and E[H] = oo, then 

^ E[I7»)| 

E[Kf E[17(2)] ■ 


(1.14) 


We think this result could be useful in other settings, particularly in other estimates involving 
local times of exclusion processes with finitely many particles. 

The paper is organized as follows. In Section 2, we give a construction of the contact 
process with stirring that allows us to separately consider the genealogy of particles and their 
motion. The construction is also very convenient to compare the process, particle by particle, 
with the associated branching process that bounds it from above. Although this construction 
and comparison were already described with words (and somewhat vaguely) in [3] and [1], as 
far as we know this is the first time that they are given explicitly. In Section 3, we show 
how Theorem 1.4 implies Proposition 1.2 and then use the construction of Section 2 to prove 
Theorem 1.1. In Section 4, we prove Theorem 1.4. 
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2 Construction of coupled processes 

Basic genealogical process. We start by giving a construction of a continuous-time branch¬ 
ing process which will be useful for coupling, and then comparing, processes of interest. Our 
construction will depend on the parameter A > 0. 

Let T be the tree defined as follows. The vertex set of T (which, by abuse of notation, is 
also denoted T), is {o} U (U^j^N”), where o is a distinguished element called the root. The 
edge set is 

^/(T) — {{o, i} . i S N} U {{(*1) • • • ) in)i (T; ■ ■ ■ 1 ini • ill ■ ■ ■ 1 in+l G N}. 

In case a = (ii,...,i„) and /3 = (ii,... ,in,in+i}i we say that a is the parent of /3 (and 
denote this by a = p(/3)) and /3 is the i„+i-th child of a (denoted /3 = Si„^;^(a)). The same 
terminology and notation is used in case a = o and /3 = i G N. 

Let Hf = {Bai Dalo-gT be a family of independent Poisson point processes on [0, oo) so 
that each Ba has rate A and each Dq, has rate 1. We regard each B^ or Da as a random 
discrete subset of [0, oo). 

We will define a process ('I't)t>o with state space {0,1, —1}'’^. Our terminology will be as 
follows: each vertex in T is called a particle] a particle is present if it is in state 1 and absent 
if in state 0 or —1. We start setting Tq = Ijo}- Now assume has already been defined up 
to time t > 0. Let 

t' = inf{s > t : s G Ba U Da for some a with Tt(a) = 1}. 

We set Tg = for every s G {t,t'); in particular, this means that, if 'l't(a) ^ 1 for every 
a, then t' = oo and for every s > t. In case t' < oo, dehne as follows. For the 

unique a such that t' G Ba U Dai consider the two cases: 

• if P G Dai we set 'I't/(a) = — 1 and = ^t{P) for all /3 7^ a; 

• if G Bai let ^ ba the smallest natural number for which 'I't(sn(Q;)) = 0, and set 
Tt/(s„(Q;)) = 1 and Tt/(/3) = for all /3 ^ Sn{a). 

It should be clear that, with this construction, the number of present particles at time t, 
: 'I'4(a) = 1}, is a continuous-time branching process with birth rate equal to A and 
death rate equal to 1. This consideration also shows that the above prescription defines the 
process (T^) for all t > 0, that is, no finite-time explosion occurs in the application of the 
recursive procedure. 

Particle positions. In addition to A, we now also fix > 0. 

We now take two more random objects, independent of each other and independent of Hj: 

• a family H^d = {Le : e G E(Z'^)} {E{lA) denotes the set of nearest neighbors edges of 

Each Le is a Poisson point processes with rate 1 on [0,00), and these processes are 
all independent; 

• a family M = {Ma{t) : a G T\{o}, t > 0} of independent random vectors of (that 
is, if («!, ti),..., (afc, tfc) are all distinct, then Ma^iti),..., Ma^.{tk) are independent). 
Each Ma{t) is uniformly distributed on the set of the neighbors of the origin of 

Given a realization of H'^d, we define a function p : {{x, s,t) : x G Z'^, 0 < s < t} —>■ Z'^, which 
we call a flow, as follows: for each x Glfl and s > 0, t >-)■ p{x, s,t) is the unique function that 
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is constant by parts and satisfies p(x, s^s) = x and, for all t > 5, 

p(x, S, t) -5, t ) if t ^ ^z\Z'^p{x,s,t—) ^{p{x,s,t—),z}^ 

yzt, xj 

p{x, s,t) = p{x, s, t-) if t G 

(we hence view t G as an “order” either to jump from a; to y or to jump from y to x, 

and p is the path obtained by following all the orders that are encountered). We also let 

Pn{x, s, t) = p{x, Ns, Nt). 


For each a G T, let 

T" = inf{t : d't(a) = 1}, T+ = sup{t : d't(a) = 1}. 

Then define 

V- ifO<t<T+; 

IA otherwise, 

where A denotes a “cemetery” state. Now assume that processes {Xa{t))t>o on U {A} 
have been defined for all a G Um=o these processes satisfy Xa{t) / A if and 

only if Ti(a) = 1. Fix a G In case t“ = oo, put Xci{t) = A for all t. Otherwise, 

fpjv (X („)(t-)+ M„(T-),T-,t) , iftG[T-,T+); 

X^{t) = ' (2.2) 

I A otherwise. 

This inductively defines Xa(t) for all a G T and t > 0, so that Xa(t) 7^ A if and only if 
'l't(a) = 1. We call Xa{t) the location of particle a at time t (with the understanding that 
absent particles are located in the cemetery state). The definition (2.2) thus means that when 
a particle appears, it is placed on a location obtained as a random neighbor of its parent’s 
location at the time, and then it moves according to the flow p^ until it disappears. 

Remark 2.1 Define the process 

= #{a G T : 'l'i(a) = 1, Xa{t) = x}, x G Z'’*, t > 0. 

Although we will not need it in the sequel, it is instructive to discuss its behavior at this point. 
In this process, particles occupy positions in it is possible that any number of particles 
occupy a single position. Each particle disappears with rate 1 and gives birth at a randomly 
chosen neighboring position with rate A (births are not forbidden at occupied sites). Moreover, 
edges of Z'’* contain jump instructions which are Poisson(N) clocks; the effect of a jump 
instruction at {x, y} is that all particles from x jump to y simultaneously, and vice-versa. 
This is the process studied by Katori in [6], inspired by the “binary contact path process” of 
Griffeath ([4])- Our task now is to construct, in this same probability space, the contact process 
with stirring this will amount, quite simply, to keeping track of particle positions and 

forbidding births at occupied sites. With this construction, at any point in time, the set of 
particles present in will be a subset of the set of particles in fit- 

Contact process with stirring. We will now construct the contact process with stirring on 
Z'^ (denoted {^t)t> 0 : with state space {0,1}^'^) and its underlying genealogical process {{r.t)t>o 
on {0,1, —I}’'*'). The construction of these processes will depend on {Hj, H^d, M) that have 
been used above and on the paths {(-^a(i))t>o : a G T} that have been defined (hence A and 
N are fixed throughout). 
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We will give a recursive definition of H that mimics that of 'h and will set 

^t(x) = #{a e T : Et{a) = 1, Xa{t) = x}, x S t > 0. (2.3) 

The construction will guarantee that, for all s > 0, 

^s(x) e {0,1} for all s > 0, x € (2.4) 

{a : S,(a) = 1} C {a : ^,(a) = 1}; (2.5) 

if ^^(a) = 1, then inf{n : Es{sn{a)) = 0} = inf{n : 4's(s„(a)) = 0}. (2.6) 


We set Ho = l{o}- Assume that (Hs)o<s<t has already been defined up to time t > 0 (and 
hence, by (2.3), {^s)o<s<t is also defined) and that they satisfy (2.5), (2.6) and (2.4) for all 
s <t. We let 

t' = inf{s > t : s £ Ba U for some a with Et{a) = 1 }. 

For s € (t, t'), we let H^ = H^ (note that this, together with (2.3), implies that ^ is then defined 
in [0, t')). If t' < oo, for the unique a for which i' S U Da-, we consider the cases: 

• if t' G Da, set Ht/(a) = — 1 and Ht/(/3) = Et{j3) for all (3 a; 

• if t' £ Ba, let n be the smallest natural number for which Ht(s„(Q;)) = 0 , and let 
/3 = Sn{a.). We consider two further cases: 

— if ^t'-{Xait'—) + Mp{t'—)) = 0 (that is, if the position where /3 is supposed to 

appear is empty), set Ht/(/3) = 1 and = Ht( 7 ) for all 7 7 ^ /3; 

— ii -{Xa{t'—)) ^ 0 (that is, if that position is occupied), set Ht/(/3) = —1 
and Ht/( 7 ) = - 4 ( 7 ) for all 7 ^/ 3 . 

These rules define H and ^ for all times. Each step in the induction preserves (2.5), (2.6) and 
(2.4), so they are satisfied for all times. Also note that 

if Hs(a) = 1 for some s, then inf{t : Ht(a) = 1 } = inf{t : 'I'f(Q;) = 1 } = T“ and 

sup{t : Et{a) = 1 } = sup{t : Tt(a) = 1 } = t+, 

that is, if a particle a ever becomes present in H, then the set of times for which it is present 
is the same for H and T. Also, 

if, for a,/3 G T, T" < < 00 and Xa(T^-) = Ap(^)(T^-) + M^(t^-), 

then Ht(/3) 7 ^ 1 for all t > 0, 

that is, a particle /? cannot appear at position x G Z"^ and time t if x is occupied at t by a 
particle a that appeared earlier than t. 

For t > 0, we denote by the smallest sigma-algebra under which the point processes 
Hf n [0,t], H-^d n [0, t] and the processes {{Ma{s))o<s<t : a G T} are measurable. Note that 
(4't,Hi,^i)t>o is then adapted to {^t)t>o- 

3 Proof of main result 

Given a G T and t > 0 so that Et{a) = 1, let n be the smallest integer so that Ht(s„(a)) = 0 
(or equivalently, by (2.6), so that Tf(s„(a)) = 0). Then consider the subgraph of T given by 

1r(Q:,t) = {a} U {j3 : the unique path in T from /3 to o intersects ( 5 ^( 0 :) : m > n}} 
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(together with all edges that are incident to two vertices in the above set). We call T(a, t) the 
set of fresh descendants of a at time t. It is a connected subtree of T. We also define, for a 
and t such that St(a) = 1 and t' > t, 

Af-{a,t,t') = #{/3 e f{a,t) : Et>il3) = 1}, J\f'^{a,t,t') = #{P G f{a,t) : = 1}. 


We observe that, for all a,t with Hi(a) = 1, 

Af^{a,t,t') < J\f^{a,t,t') 


(3.1) 


and, since [Af'^{a,t,t + s))^^^ is a branching process with birth rate A and death rate 1, we 
have 

on {Ht(a) = 1}, for all t' >t, K {a,t,t') \ 5'f] = (3.2) 

We now define two classes of events, I{a,t) and J{a,t), which will both correspond to 
situations where there is a 0 —)• 1 transition in which does not occur in H. We will be able 
to guarantee that, if one of these events occurs, the number of particles present in S is strictly 
smaller than that in 'I'. Both definitions will depend on the following quantity: 


log 


(3.3) 


We start with the definition of I{a,t). Fix a G T, t > 0 and let n be the smallest integer 
such that St(sn(Q;)) = 0. Let /3 = Sn{a) and 7 = Sn+i(a). Assume that the following occurs: 

Et{a) = 1, {Da U D/sU D^U BfjU B^) r\[t,t + t*] = 0 , if{Ba r\[t,t + t*]) = 2. (3.4) 

Additionally, letting Ba n [t, t + t*] = {Ti, T2} with Ti < T2, assume that 

Xa{T 2 -) + M^{T 2 -) = Xp{T 2 -). (3.5) 

(This implies that either j5 obstructs the appearance of 7, or some earlier particle obstructs 
the appearance of both /3 and 7). Let I{a,t) be the event described in (3.4) and (3.5). We 
then have 

on I {a, t), Af^{a, t,t + t*) < 2 and A^'^(a:, t,t + t*) = 3. (3.6) 

We now turn to the similar dehnition of J{a,t). Again hx a G T, t > 0, n the smallest 
integer such that Et{sn{a)) = 0 and /3 = Sn(a). Also let 7' = si(/3). Assume that 

Et{a) = l, {Da 0 D 0 [JDY 0 By)n[t,t + t*] = 0 , 

if{Ban[t,t + t*]) = #{Bf^n[t,t + t*]) = l. 

Letting Ba n [t, t + 1*] = {Ti} and Bp n [t, t + t*] = {T 2 }, we will also require that 

T,<T 2 , Xa{T 2 -)=Xp{T 2 -) + My{T 2 -). (3.8) 

(This implies that either /3 obstructs the appearance of 7^, or some earlier particle obstructs the 
appearance of both j3 and 7'). We define J(a, t) as the event specified by all the requirements 
in (3.7) and (3.8), and note that 

on J(a, t), Af^{a, t,t + t*) < 2 and {a, t,t + t*) = 3. (3.9) 


Due to the invariance under time shift of Hj, and M, 

on {St(a) = 1}, P[/(a,t) I = P[/(o,0)] and P [J(a, t) | 5i] = P[J(o, 0)]. (3.10) 
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Also, inspecting the above definitions one can show that 

P[/(o,0)] =P[J(o,0)]. 

We are now ready for the key estimate of this section. For any t > 0, 


(3.11) 


E[#{aeT:Ht+t*(a) = l} I = E 


a€T:Et{a)=l 


dt 


(3.1),(3.6),(3.9) 

< E 

Y2 Af^ {a,t,t + t*) 

dt 

-E 


dt 


oeT:St(o)=l 



aGT:St(a)=l 



( 3 . 2 )^ 3 . 10 ) ^ ^ ^ ^ -P[/(o,0)] -P[J(o,0)]) 

) #{a G T : Et{a) = 1} • - 2P[/(o,0)]). 

Applying this estimate recursively, for any /c G N we have 


E 


#{x G : 4t*(x) = 1} ^ = ^ E[#{a G T : = 1}] < - 2P[/(o,0)])' 


This shows that 

if e‘‘A-i) -2P[/(o,0)] < 1, then 

P[^i / 0 for all t] = lim P[4i* 7 ^ 0] < hm E[#{a; : ^kt*{x) = 1}] = 0. 

fc—>-oo k^oo 

We now estimate P[/(o,0)]. We denote (3 = si(o), 7 = 52 ( 0 ) and 

E = {{Do U Dp\J D^U B^) n [0,^] = 0, #(^0 0 [0,t*]) = 2}, 

so that ^ 

P[^] = l^^exp{-3r - 3Xt*}. 


(3.12) 


(3.13) 


On E, let Ti < T 2 denote the two elements of Bq H [0,t*]. Let / denote the joint density 
function of Ti,T 2 conditioned on E; then, f{ti,t 2 ) = ■ 1 jo< 7 <i 2 <t*}. Now, 

P[/(o,0) \E]= nMfi{Ti) = zi, M^(T2) = Z 2 ] 


Z\,Z2&3-. 

zi,22~0 


t* rt* 


S.i I- 


(t *)2 


P[pw(0,0,ti) = X, pN{0,0,t2) = y, PNix + Zi,ti,t 2 ) = y + Z 2 ]dt 2 dti 


Y] P[M^(7i) = Zi] • [ [ ■F[pN{0,0,t2-ti) ^ pN{zi,0,t2-ti)]dt2dti 

Jo Jti v* ) 


zi'^0 


L.r r 

2d Jq Ju 


2d Jo Jh {t*V 


P[pAr(0, 0, t 2 - h) ~ pNiz, 0, t 2 - tl)] dt 2 dti 
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for any fixed z ~ 0 in Changing variables twice, this is further equal to 
^ IP 0, s) ~ pn{z, 0, s)] • it* - s) ds 

1 / u \ 

^ dFN Jo ^ “ ¥n) 

We now fix e S (0,1) whose value will be chosen later. The expression in (3.14) is larger than 

2 fet*N 

^^•( 1 -e) P[pi( 0,0 ,m) ~ pi(z,0,'u)] d-u. 

= • P [|{m < et*N : pi{0, 0, u) ~ pi{z, 0, u)}|] . (3.15) 

We abbreviate 

g{t) = E[|{m < i ; pi(0, 0, u) ~ pi{z, 0, n)}|]; 

by Proposition 1.2, to be proved in the next subsection, we have limt_;.oo ^ 

d > 2. Putting together this fact, (1.8), (1.9), (3.13) and (3.15) we see that, if N is large 
enough (depending on e), 


P[/(o,0)] > • g{N) ■ • exp{-3C - 3Xt*}. 


If we further assume that A = Aw 


dt*N 

N^oo 


> 1 and use (3.3), then for N large enough we have 
3 t*g{N) 


P[/(o,0)]>(l-e) 


2dN 


(3.16) 


We are now ready to conclude. Assume that A = 

N^OO 


particular, by by Proposition 1.2, we have Aw 


= Aw = 1 + for some 0 < 1 (so in 

> 1). Then, if N is large enough. 


gt*(A-l) _ < (1 +e) • 




t*Q(N) (3.16) 

-%r < 2P[/(o,o)], 


where the inequality (*) holds if e is small enough that > 6. The desired result now 

follows from (3.12). 


3.1 Exclusion dynamics of two particles: proof of Proposition 1.2 

We will denote the vectors in the canonical basis of Z'^ by ei,..., e^. We write 

M) = {ei, —ei, ■ ■ ■, Crf, —Gd}, ATo = A/q U {0}. 

For d > 3, as in the Introduction, we denote by G{x, y) the Green function of discrete-time, 
simple random walk on Z'^ and by {Ag, Bs)s>o the process given by the positions of two 
particles moving on Z'^ under exclusion dynamics with rate 1 per edge of Z"^. This means that 
{As, Bs)s>o is the Markov process on (Z“*)^ with generator 

(£exc/)(a, b) ='P2 [^{a+z^b} ■ if {a + z,b)- f{a, b)) + t{b+z^a} ■ if {a, b + z) - f{a, b)) 
zSA/'o 

~^^{a+z=b}U{b+z=a} ' (/(^) ®) /(®)^))] ■ 
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We always assume that Aq — Bq G A/q, that is, the particles are initially at neighboring 
positions. 

Now let 

Xs = As — Bg, s > 0; 

observe that Xg G Z‘^\{0} for all s and that Xq G A/q. The generator of {Xg) is 


, |2EzeAro(/(® + ^) ifx^Afo; 

[Cxf)[x) = < 

[ ifi-x) - fix)) + 2 Y)zeMo\{-x}ifix + 2 ) - fix)) otherwise. 

For comparison purposes, it will be useful to also consider the process iYg)g>o, which has 
generator 

(£y/)(x) = 2 ^ ifix + z)- fix)); 

z£j\fo 

we also assume Yq G Aq. (Yg) is hence a continuous-time simple random walk on that 
jumps from each vertex with total rate dd. (Al^) and (Yg) have the same behavior except 
when iXg) is in some position z G A/q, in which case, instead of jumping to 0 with rate 2, it 
jumps to —z with rate 1. 

Now let 

= inf{t >0:Xt(f A/'o}; 

= inf{t >0-.Yt^ ATo}; 

UYXo = e [0, : Yt G Mo}; = |{t G [0, U^] : Yt = 0}, 

so that U'^ = I/TA/q _j_ gy ^ simple computation involving geometrically distributed 

random variables, it is easy to check (or see the proof of Lemma 3.4 in [1]) that 

E = E ^ 

Proof of Proposition 1.2. First assume d > 3. Let q denote the probability that iXg) never 
returns to A/'o after first leaving it. Note that this is the same as the probability that iYg) 
never returns to A/’o after first leaving it. We then have 


E[|{s<oo:X, gA/'o}|] 

_ K[U^] (3U7) E[U^Mo] 

Q Q 


E[|{s<oo:y, gAAoII] 


4 G(0,0 = ^(0(0,0)-1), 

zeXo 


where the last equality follows from G(0,0) = 1+^ EzeA/b 0) ~ 1+G(0, ei), by symmetry 
of the Green function. 

For d = 2 we have 


E[|{s<t:X, gAAqII] ^ E[|{s<t:y, gA/'oII] _ E [|{s < t : X, G AAo}|] ^ 

logt logt E [|{s < t : Y; G AToII] ’ 

Now, Theorem 1.4 implies that 


y E [|{s <t:Xse A/’o}!] _ E [U^] (3U7) 2d _ 4 
E [|{s < t ; y, G :Eo}|] ~ E[UY] ~ 2d + l ~ 5' 


(3.19) 
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Now let us treat the first quotient in (3.18). By the Local Central Limit Theorem (Theorem 
2.5.6 in [8]), for any 2 : G Z'^, 

= Z] S—J-CXD .. 

where ft(x) = is the probability density function of a Gaussian vector {y,W) so 

that V and W are independent and have mean 0 and variance t. Hence, 


lE[|{s<t:y, g:Fo}|] 


logt 


1 /"* TO) rw 1 u 2 d + 1 

- - > / P Ts = -z ds-^- 

logt ^ io ^ ^ 

zeAfo 


27r • 4 87r 


= — (3.20) 


and similarly, 


P [Ks < t : G A/qII] t—>-00^ 2 d 


1 


log t 27r • 4 27r 

Now, (3.18), (3.19), (3.20) and (3.21) imply the desired result. 


(3.21) 


4 Two-type renewal processes: proof of Theorem 1.4 

Take (K)n>o, {Un^)n>o, {Sn^)n>0 and for i = 1,2, as in (1.11), (1.12), (1.13). Theorem 
1.4 follows immediately from the two lemmas: 


Lemma 4.1 = 1, = 1 and E[(t/*^^i)^] < 00 , then 

hm -= 1. 

E[k\ >] 

Lemma 4.2 Let > 0. and E[y] = 00 , then 

E[^ ^ ^ 

E[Kf ' 


(4.1) 


(4.2) 


Before we prove these results, let us give a definition. For i G {1,2} and t > 0, we let be 

(i) (i) 

the unique value of n such that <t< 82 ^+ 2 '^ note that 


7V«-1 

n=0 


t, s 


(0 

2 N, 


(i) 



5 


(0 

2 N, 


(i) 


= t 


7V«-1 

Vn — max 

n=0 



t-S 


(0 

2 N, 


(i) 



(4.3) 


(4.4) 


(Here and 
zero). We 
we have 


in what follows, we interpret sums of the form Yln=a ''^ifh a > 6 to be equal to 
also observe that, since and V are assumed to be positive random variables, 


hm = 00 and 

t^OO 

a.s. 


lim E[A'"j^*^] = 00 . 

t^oo 


(4.5) 


We also claim that 

if E[t/(*) + H] = 00 , then lim -E[N'i^*^] = 0. 

t—>-oo t 


(4.6) 
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Indeed, the fact that Nt/t —)• 0 almost surely as t —)• oo follows from the Law of Large Numbers, 
so it is sufficient to prove that supj>]^ E[(Nt/t)^] < oo. This can be shown by fixing e > 0, 
(5 > 0 such that + V > S\> e and noting that 


< min ^ m > 0 : >«} > rt/41 \; 


m—1 


n=0 


the right-hand side is distributed as a sum of /5] geometric random variables with parameter 
larger than e; this gives < |~|] • Hence (4.6) is proved. 

Proof of Lemma f.l. Define 


m—1 


m—1 


=max It - t'l : ^ [/(D <t<Y,U^n\Y. ^ ^ E 


m> 0. 


y n=0 n=0 n=0 

We will now prove two claims. 

Claim 4.3 There exists C > 0 such that, for all k > 1, 


n=0 


E 


max Am 

0<m<k 


< Ck-i. 


(4.7) 


To prove this, start noting that, since = 1, Am is equal to 


max • 


m — 




n=0 


m—1 


m 


+ 1 - E 


n=0 


< max • 


m 




n=0 


m—1 


m — 


T.vi? 


n=0 


+ 1) 


SO that 


max Am < max 

0<m<fc l<m</c+l 


m—1 


m — 




71=0 


+ 1 . 


( 2 ) 

Since the Un ' are independent and identically distributed with mean 1, the Reflection Prin¬ 
ciple ([8], Proposition 1.6.2) implies that 

max]^<m<fc+i A m 


E 


fc3/4 


< 2-E 


fc + i-EtoL'’ 






P/4' 


Write Wk = 


k + l- Z)n=0 


( 2 ) 


Wk 

hmsup^=== 
k^oo Vkloglogk 

a.s. 


. By the law of the iterated logarithm, 
< oo hence 


Wk 

lim twtt = 0 . 


fc—foo 


(4.8) 


Then, 


E 


Wk 

*3/4 


< E 


Wk 


fc3/4 ■ 


+ E 


(my 

V*3/4y 


Wk 

= E -A . 1 






Var(Wfc) 

*3/2 


= E 




k + l 


The expectation on the right-hand side vanishes as * —)• oo by (4.8) and the dominated 
convergence theorem. This completes the proof of the claim. 
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Claim 4.4 


(2) ^(1) 


Kl — K- 


< max 

n -m A/" ^ ^ ^ 


t > 0. 


(4.9) 


To show this, fix t > 0. First consider the case then, it readily follows from (4.3) 


that 


(2) ^(1) 


k; — K 


< A^(i)_^. Next, assume that 




We then have 


/(i). 


r(2) 


hi) (Ilf) 


k; = t — 


t — Vn — max ( 0, t — 


n=0 


:(i) 


/(!)_ 


(4.10) ^ (4.4) . . 


n=0 


Hence, 


j(2). 


(2) ^(1) 


k; — K 


- .( 2 ) _ .( 1 ) ( 4 - 1 - 3 ) 

— tv^ rvj- — 


n=0 


42) 


C/,(2)+mm(^t, 

C/W - Linft, 




(4.10) ^ * , 

< Y1 

n=0 


( 2 ) 




n=0 


-s 




-s' 


2 nSSi 


( 2 ) 

AN 


( 2 ) 


7(2). 


( 2 ) 

2N, 


( 2 ) 


E 


n=0 


^ ^A7(2)_ 


(4.10) 

< max Am ■ 

n m A7 (^) 


-5 


( 1 ) 

2 N^ 


( 1 ) 


(4.10) 


Now, a symmetric argument shows that, when A^^) < Aj^\ we have 


(2) ^(1) 


k; — K- 


<A.,w_ 


(4.10) 

< max Am ■ 


This completes the proof of (4.9). 

We are now ready to conclude. We have 


e[ac: 


(2)i 


EK 


( 1 )- 


-1 


= E[4')]-1- E[icS")]-E[kW] <E[.cS')]-4-E 


max A^ 
mW 


(4.11) 


Now note that A^^^ only depends on (Fn)n>o and (A„)„>o only depends on {Un^)n>o- Hence, 
A^^^ is independent of (A„)„>o. Thus, 


E 


max Ar 


0<m<N^ 


( 1 ) 


E*^ 

k=0 


max Am I a/^^ = k 

0<m<k 


Af ^ = k 


E"^ 

k=0 


max Am 

•P 

Af ^ = k 

0<m<k 


- 


< C^A:4 .p 


k=0 


Af ^ = k 


( 1 ) 


= CE 




(Hit 


Since A^^^^ —)• oo almost surely as t oo, we have E (A^*'^^)4 /E A^ 
shows that the right-hand side of (4.11) vanishes as t ^ oo, completing the proof. 


0 as t —)• oo. This 
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Proof of Lemma 4■2- If the result follows from Lemma 4.1 and a change of scale. 

So we assume (without loss of generality) that It follows from the assumption 

that and (4.3) that 

■Nf'^\ < a^®), i G {1,2}, t > 0, 


so, writing 


E[Kfa(2)E[ivf+ E[Kf'> - a(2)ivf)] 
and using (4.5), we see that (4.2) will follow once we prove that 


mt 


(i)i 


lim 

t^oc E[ivf)] 


= 1 . 


(4.12) 


n=0 


With this in mind, we now proceed to bound Observe that 

k-l 

< k if and only if t < Vn- 

Since 

„(2)jvf) + ^ ^ t}, = s^^a) < t, 

n=0 n=0 

applying (4.13) with i = 2 and k = implies that We now claim that 

A/'f ^ - A^f ^ < 1 + inf < 


(4.13) 








m 


> 0 : • m + ^ 14 > (a^^^ — «^^^)(A^i+ 1 ) 


n=Ar)®4l 


> . (4.14) 


Indeed, assume m belongs to the set of which the infimum is taken on the right-hand side. 
We have 




N, 


( 1 ) 


/ 


a(2 )(ivW+m + l)+ ^ Vn = a^'^'>Ni^'^+ Y,yn + oi^‘^'^ + 


n=0 


n=0 




\ 


Qp‘'^m 


\ 


+ E 

n=7v)®4l / 


> a(2)^a) + ^ v; + + (a(^) - a(2))(Ar^(i) + i) 


n=0 


N 


( 1 ) 


= aW(ArW + l) += 


27v)®42 


n=0 


> L 


so, applying (4.13) with i = 2 and k = -|- m -|- 1 we get -|- m -|- 1, proving 

(4.14). 

Now, if Vq, V(, ... are random variables distributed as V and independent of {Un^)n>o and 
(I4)n>0! fl^e right-hand side of (4.14) has the same distribution as 


m—1 


1 -|- inf < m > 0 : -|- J]yXa(i)-a( 2))(jvf + l) 


n=0 
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Conditioning on (and denoting by the distribution of (4.14) then yields 


roo 

E[|ivf-iv(')|]< / E 

Jo 

roo 

= / E 
Jo 


m—l 


1 + inf < m > 0 : ^ > (a^^^ — a^^^)(s + 1) 


1 + iV, 


( 2 ) 


n=0 

fJ'Mwids). 


(ds) 


Now hx e > 0. By (4.6), there exists C > 0 such that E[A^i^^] < C + es for all s > 0, so the 
above gives 

E[|ivf) < 1+ C + e-E[(a(i) - a(2))(iV^(^) + 1)]. 

Finally, by (4.5), 

0 < limsup -TTT^-^ < e(a^^ - 0; 

i^oo K[NI ’] 

since £ is arbitrary, this completes the proof of (4.12). I 
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